Interaction engineering for environmental probing 
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We study the conditions for the probing of an environment affecting one party of a bipartite system of inter- 
acting objects by measurements operated only on the other element. We show that entanglement plays no crucial 
role in such an environment-characterization. On the other hand, if an interaction is established between the two 
parties, information can be reliably gathered. This result holds for both discrete and continuous variables and 
helps in the interpretation of recent experiments addressing the properties of mesoscopic objects. 
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"How much of the dynamics of an open system can be re- 
vealed without directly accessing it?" This question is becom- 
ing increasingly relevant to the performance of experimental 
quantum information processing due to the development of a 
new generation of protocols dealing with systems which are 
only partially accessible from outside. Built-in charge qubits 
interacting with the field of on-chip superconducting waveg- 
uides 1 1], electron spins in quantum dots J2], micro- or nano- 
electromechanical systems and their optomechanical counter- 
parts and single-molecule junctions [5] are striking ex- 
amples of the situation depicted above, where only partial ac- 
cess to the components of a register is possible. Another sce- 
nario where a similar question holds is given by many-body 
systems of (possibly interacting) particles, such as optical lat- 
tices [6], where the direct addressing of a specific element of 
the system is difficult because of too-close inter-element dis- 
tances. In this case, a winning strategy for the revelation of 
the dynamics of one of the particles in the register (regardless 
of its nature) can be the use of controllable ancillary devices 
acting as detection stages. The inability of directly addressing 
a system is then bypassed by coupling it to an ancilla (not nec- 
essarily being a qubit) that is easier to manipulate and detect. 

Our question has natural pragmatic relevance and is cer- 
tainly experimentally stimulating. We show that it is also 
dressed of fundamental importance, as it can be related to the 
role played by entanglement in the revelation of the dynamics 
of a system. By studying discrete and continuous variables 
(CV's), we show that the engineering of "interactions" be- 
tween the detecting stage and the system to study, rather than 
entanglement, is crucial in an environment-revelation process. 
This holds under general conditions: Entanglement alone is 
not sufficient for indirect environment revelation. Even in the 
presence of entanglement between (non-interacting) system 
and detecting stage, tracing the system before or after its in- 
teraction with the environment is irrelevant, when looking at 
the detecting part alone. In our analysis we use a minimum- 
control approach, where no fast/precise interaction- switching 
is required and the environmental inference is performed sim- 
ply by looking at the equilibrium state of the system, which 
is in general easily measured. These features allow for an 
original reinterpretation of recent experiments performed on 
mesoscopic optomechanical devices (4j]. 



Qubit case.- In order to tackle the central question of our 
study, we start considering the simple case of two qubits, la- 
belled 1 and 2 and characterized by the respective transition 
frequency ujj (j = 1,2). The free evolution of the qubits 
is therefore ruled by the Hamiltonian Hf = uj\a z ^/2 + 
^2^2,2/2 (we assume h = 1) with the z-Pauli matrices 
°z,j = - and {l°> A 1 )} 3 the computational 

states of qubit j. Qubit 2 embodies the system whose dy- 
namics we want to characterize through the detecting stage 
represented by qubit 1. In order to account for the case where 
qubit 2 is originally part of a multipartite register whose state 
we are not interested in, we allow for 2 to be prepared in a 
general mixed state. This would be the case, for instance, 
when this qubit is originally entangled with the rest of a reg- 
ister. The second assumption we make is that qubit 1 is in 
a pure state, well isolated from the environment affecting 2. 
Detection- windows can be opened on 1, during the dynamical 
evolution. This assumption does not affects the generality of 
our results and can be easily relaxed. 

We consider a non-switchable interaction between 1 and 2 
described by the coupling Hamiltonian Hi. As qubit 2 and 
its surrounding physical system are unavailable to direct in- 
vestigation, we assume no knowledge about the environment 
B affecting qubit 2. The only working condition we take is 
that B is Markovian. This simplification, allowing us to treat 
the problem with a Liouvillian approach, holds in many phys- 
ical scenarios such as high-temperature circuit-QED, where 
the Bloch-Redfield formalism must be pushed beyond the 
standard secular approximation and "quantum optics" mas- 
ter equations (ME's) become the most appropriate tool 10]. 
When this assumption cannot be made (as for the case of 
long-memory nuclear spin baths in semiconductor quantum 
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FIG. 1: Sketch of the physical situation. System 1 is completely 
accessible while non-accessible system 2 is embedded in an environ- 
ment, 23, whose parameters T, 7 are unknown. Systems 1 and 2 can 
interact via a non-switchable coupling Hi. 
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dots), existing tools such as effective Hamiltonian approaches 
or non-Markovian ME could be used @]. 

The choice of Hi is a setup-dependent issue. Its form is 
usually guided by naturally (or easily) realized interactions 
specific of a chosen implementation. However, our scheme 
is explicitly hybrid and thought as the result of the combi- 
nation of register and detection stage having different phys- 
ical nature. As such, it is reasonable to assume the ability 
of engineering the most suitable form of coupling, possibly 
by changing the type of detection stage being used. To illus- 
trate the main point of our analysis, we find instrumental the 
study of two particular models. The first is an anisotropic xy 
coupling, which has been largely used in problems of quan- 
tum dynamics in spin chains. We take it because it natu- 
rally arises in quantum electrodynamics-like systems where 
two qubits dispersively interact with a bus via the off-resonant 
Jaynes-Cummings model (9|]. This latter is almost ubiqui- 
tous where a spin-boson interaction is arranged 1101 . The 
xy model reads H xy = J x Vx,i&x,2 + Jy&y, i&y,2- The sec- 
ond model we consider, H c bf = (#/2)(H2x2 + ^,1)^,2 » 
describes a bit-flip of qubit 2 conditioned on the state of 
qubit 1. Its use is related to the extension to CV's that 
we perform in the last part of this paper. Here, J XyV and 
g denote coupling strengths. The non-unitary part of the 
evolution is entirely ascribed to qubit 2 being exposed to 
its environment. We consider a general case of B being 
both dissipative and dephasing. The Liouvillian Cb(q) act- 
ing on the density matrix g of the system is therefore 
C b (q) = -r(n+l)({|l) 2 (l|, e }-2|0) 2 (l|e|l> 2 <0|)- 
rn({|0) 2 (0| ,q}-2 |1) 2 (0| g|0> 2 <1|) - 7[**,2,**,20]. The 
first two terms account for a dissipation occurring at a rate 
T (n is the average thermal occupation number of the envi- 
ronment), proportional to its temperature, while the last term 
describes dephasing at rate 7. With this notation, the evo- 
lution of the density matrix of the system is given by the ME 
d T g = -i[H t , q]+Cb(q) with 7^ = H f + Hi andHi = H xy 
(Hi = Hcbf) for the choice of the first (second) model. H xy 
can be written as H xy = <j +5 i(J<j_ 5 2 + ^+,2) + h.c, where 
J = J x + J y , S = J x - J y and a±j = (a x j ± ia y j)/2. 
For 5 7^ 0, [Hxy, J2j &z,j] 7^ that prevents the total num- 
ber of excitation to be a good quantum number. If we con- 
sider the cbf model, instead, we have [H c bf,&z,i] = 0, im- 
plying that the Hilbert space is partitioned in two sectors, 
each spanned by states having qubit 1 in its ground or ex- 
cites state, respectively. The eigenstates of H c bf are |00) 12 
and |01) 12 for the first sector (eigenvalues — ft/ 2 and —u>/2) 
and [-(u 2 /2g) =F y /l + (^/2g ) 2 } |10) + |11) (with eigen- 
values (oji =p 002 yl + 4^/cJf)/2) for the second. Here, 
£1 = u)\ -\- CO2 and O = oui — 002- Obviously, to get entan- 
glement through H c bf, we should start with a superposition 
costflO)! +e i *sintf|l) 1 (& e [0,tt/2], <p e [0, tt]). 

Independently of the interaction being considered, the dy- 
namics of the system is best evaluated by projecting the ME 
onto the elements of the two-qubit computational basis, i.e. 
by taking d t Qi m , P q = -i (lm\ [H u q] \pq) + (lm\ t B (o) \pq), 



where Qim, Pq — (Itti\ P \pq) (l,m,p,q = 0, 1) is an element 
of the density matrix of the two-qubit system. The explicit 
form of this set of coupled differential equations is straight- 
forwardly found by direct calculation. 

By introducing the vector of density matrix elements Vi 
(i = xy, cbf) such that (vi) p = Qp with p = 0, .., 15 and 
p the same number expressed in binary notation, it is possi- 
ble to write such B loch-like equations in the vectorial form 
d t Vi = E^Vi, where 5^isal6xl6 matrix of coefficients 
that is easily extracted from the Bloch equations. The dynam- 
ical evolution can be cast into the form Vi(t) = e Si£ Vi(0), 
where the exponential matrix includes the effects of B. It is 
straightforward to see that, in general, the steady state of the 
system, independently of the coupling model being assumed, 
can be determined as the eigenvector v| s corresponding to the 
null eigenvalue of 5V For the xy coupling, this immediately 
leads to q%% = (T/d xy )g s x s y + (J 2 5 2 GN / 'd xy )t 4x4 . Here 
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with N = 2n+ 1, G = TN + 27 and, explicitly, A± = (N =F 
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l )S 2 + (N ± 1)(G 2 + n 2 )J 2 , S = -iTJNS 2 (G - 



iO), T = iTJ 2 S(G 

uo 2 )5 2 + (G 2 + Q 2 )J 2 } + 4GJ 2 5' 2 }. The steady state g s xy 
thus a mixture of Bell states, its separability resulting from 
the trade-off between the parameters entering such a mix- 
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d~ h f = TrB. It is easy to see that is separable. A way 
to understand it is the following: If we choose a state having 
$ = 0, for instance, we know that the state has to remain sep- 
arable throughout the entire dynamics leading to the steady 
state. However, the properties of the steady state do not de- 
pend, by definition, on the specific instance of initial state be- 
ing considered. Therefore, the only possibility we have is that 
the steady state is always separable. Of course, dynamically, 
there can be a transient during which quantum correlations are 
present. This is a point of interest of the following analysis. 

The entanglement measure that we adopt here is based on 
negativity of the partial transposition criterion, which is nec- 
essary and sufficient in the general two-qubit case 111211 . The 
dynamical entanglement Ni(t), calculated via the knowledge 
of Vi(t), can be compared with entanglement at the steady 
state Nf s . Figs. [2] show some significant examples of such a 
comparison. In panel (a), relative to the xy model, N xy (t) is 
persistent in time, stabilizing at about 0.06. This is in con- 
trast with the situation of panel (c), where N c bf(t) (for qubit 
1 prepared in (|0) 1 + |l) 1 )/V / 2) goes to zero. Obviously, in a 
"warm" environment having n = 1, N xy (t) is no more persis- 
tent in time. This is better shown in Fig.[2](b), where N xy is 
plotted against n and shown to be almost linearly decreasing 
(see the figure caption for more details). As anticipated, 
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FIG. 2: (a): N xy against Jr (solid line) for J = 0.3, S = 0.1, u) = 
0.3, Q = 3,7 = 0.001, T = 0.1 andn = 10" 3 . Qubit 1 is prepared 
in |0) 1 . We took 100 initial density matrices for qubit 2, each cor- 
responding to a Bloch vector extracted from a Gaussian sample. We 
plot the average entanglement of such an ensemble (solid line). The 
dashed line shows N*y. (b): N* y against n for the same parameters 
in (a), (c): N cbf for g = 0.3, uj = 0.3,7 = 10" 2 , r = 1, Q = 3 
andn = 10" 3 . Qubit 1 is prepared in (JO^ + |l> 1 )/>/2. 



is always zero. For the cbf Hamiltonian, in Fig.[2](c) we show 
an example addressing a strong environment. An emerging 
feature is that the time at which the system reaches its steady 
state depends very weakly on the initial state of qubit 2 but 
critically on the parameters of its environment. The larger 
(7, T, n) with respect to (J, S) or g, the quicker p| s is set, this 
being entangled depending on the trade-off between coherent 
and incoherent dynamics. 

Environment probing.- Motivated by recent efforts pro- 
duced, for instance, in circuit-QED |yj], in order to discern the 
effect of the coupling between qubit 2 and its environment, we 
consider the spectrum of emission from qubit 1. The choice 
of this figure of merit, which is standard in quantum optics, 
is strengthened by the fact that "easy" single-qubit measure- 
ments are required to reconstruct it. We aim at determin- 
ing the spectrum S(u) = 3? [J °° dre^ r (a +) i(r)(j_ ) i(0)) ss ] , 
where (<j +5 i(t)<t_ 5 i(0)) ss is a two-time correlation function 
evaluated in the steady state. The linearity of the B loch-like 
equations allows the use of the quantum regression theorem 
to obtain the two-time correlation function we need [11]. We 
indicate with C^(r) the vector of two-time correlations rel- 
ative to a choice of Hi. Through the corresponding Bloch- 
like equations, its time-evolution is ruled by d T Ci(r) = 
M.i Ci(r) with M$ a square matrix of coeffcients and C(0)$ 
easily determined via gf s . By arranging each C%{r) so that 
(Ci(r))i = (<3-+,i(t)<j_ 5 i(0)) ss , we arrive at the elegant ex- 
pression S(v) = (3? [(ivt - M)- 1 C(0)]) 1 . Analytic solu- 
tions can be gathered and turn out to be identical to what is 
found by brute force solution of the dynamics of C^(r)'s and 
Cramer's rule HH [3. It immediately appears that, regard- 
less of the interaction model being adopted, a discrimination 
of the parameters characterizing the bath B affecting qubit 
2 is possible. Indeed, let us examine some particularly well- 
visible situations. Fig. [3] (a) addresses the xy model, showing 
the case of two frequency-degenerate qubit in the presence of 
a weak bath affecting qubit 2. Both dissipation and dephasing 
are considered. Two Lorentzian peaks appear at small tem- 
peratures of B, with a dip in correspondence of the common 
frequency uo\ =002- This is a clear effect of quantum interfer- 
ence, still possible in such a quasi-unitary situation. As soon 
as the thermal nature of B increases, the incoherent dynamics 
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FIG. 3: (a): S(v) against frequency v and bath temperature n for 
the xy model and J 0.3, 5 0.1, uo = 0, 7 = F = 10" 3 , Q 
3. (b): Effect of dephasing on the emission spectrum for the same 
parameters used in panel (a) and n — 10 -2 . 



reduces the dip and merges the Lorentzian peaks into a sin- 
gle thermal peak. The amplitude of the two Lorentzian peaks 
depends on T while the dephasing rate 7 has strong effects 
on the peak-merging process. For a very small temperature, 
an increase in 7 makes the interference effect quickly disap- 
pear (see Fig. [3] (b)) due to the cancellation of off-diagonal 
elements in g s x s y . On the other hand, larger dissipation implies 
that the amplitude of the peaks fades more rapidly (an increase 
in T from 0.001 to 0.01 needs a reduction of the n range from 
300 to 30). That is, the engineered interaction enables the 
discrimination of different environmental conditions to which 
qubit 2 is exposed. We can infer whether B is dissipative (de- 
phasing) only: for instance, 5(1.5) at T = 0.15, n = 10 and 
7 = 10 -3 is roughly six times smaller than the value corre- 
sponding to the case of T <-> 7. We can have even more insight 
and discern if any asymmetry is present in the free dynamics 
of the qubits. Indeed, for Co ^ 0, the spectrum becomes asym- 
metric, with one of the peaks more pronounced and fading 
more slowly than the other one. Analogous features are found 
in Figs.Hl where the cbf model is studied. 

In the case of Fig. 0(a), N X y disappears for n > 0.03 while 
for Figs. HI strikingly, we know that g s c lj is never entangled. 
Nevertheless, the process of environment discrimination is ef- 
fective. To fix the ideas, let us concentrate on the xy case: the 
g s x s y associated with Fig. 0(a) is not diagonal but keeps (small) 
coherences. Although at most classical correlations are shared 
by the qubits (as witnessed by the fact that the mutual infor- 
mation in g x s y , accounting for the "total" correlations between 
the qubits, is non-zero), as long as the interaction is present 
(and not overcome by decoherence killing even the classical 
correlations in the system), the environment can be character- 
ized. Indeed, for the xy model with J XjV = (cbf model with 
$ = 0), no inferrence from the observation of 1 is possible 
(the emission spectrum is flat), even in the presence of pre- 
constituted entanglement. However, any small deviation from 
this condition brings the probing mechanism back to effec- 
tiveness, regardless of the absence of entanglement. Only the 
interaction rules the possibility of probing the environment. 

CV case.- The results described so far can be extended 
to the CV scenario without any changes to our conclusions. 
Rather than slavishly reproduce the approach followed in the 
qubit case, here we show that entanglement does not imply 
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the ability to infer the properties of an environment affecting 
an inaccessible part of a bipartite CV system. Let us con- 
sider a two-mode squeezed state of squeezing factor r (the 
modes are labelled 1 and 2). While mode 1 undergoes a uni- 
tary free evolution, mode 2 interacts with a bath that, without 
affecting the generality of our discussions, is taken as dissi- 
pative and characterized by the (unknown) parameters (T, n). 
No interaction is assumed to connect the two modes. This 
restricts the study to a dissipative dynamics best described 
through the Markovian Fokker-Planck equation (FPE) for the 
Wigner function of the two-mode state 111411 . The solution of 
the FPE is the convolution of the initial Wigner function of 
a two-mode squeezed state Ill4ll and that of a thermal state 
modelling the environment [ 15] . In details W(a,/3, r) = 
A/"exp[-^f |a| 2 - ^-\(3\ 2 + + a*/?*)], where Af is a 

normalization factor, T = y/l — R 2 = exp[— Tr/2], D = 
2T 2 + 2R 2 N cosh 2r . Here, E a = T 2 cosh 2r + NR 2 , Ep = 
cosh2r and F = Tsinh2r while (3 (a) is the complex 
quadrature variable of mode 2 (1). By integrating W(a, (3, r) 
over /3, we obtain the thermal Wigner function Wi(a, r) = 

2| |2 

(2/irEp) exp[ — which carries no information about the 
environment affecting mode 2, confirming that entanglement 
is not sufficient for environmental probing. 

In order to address the case of two interacting CV's, we 
consider a non-trivial case of experimental interest. Let us rea- 
son a posteriori by considering the radiation-pressure Hamil- 
tonian H rp = x&l&i(&2 + ^2) responsible for the coupling 
between a field mode and the vibrating mirror of an opti- 
cal cavity (modelled as a harmonic oscillator). Here, a\ 
(0,2) is the annihilation operator of the field (mirror) and x 
is the coupling constant 

Elm Hi. 

The time-propagator 
U rp = Q-^tnaih^ +h) ^ s a displacement operator 111411 whose 
complex amplitude depends on the intensity of the field. If 
the Hilbert space of mode b is restricted to the subspace with 
at most one excitation, U rp is isomorphic to e~ %ncbft . H rp 
is thus the CV extension we are looking for. The Bloch-like 
equations we have discussed in the qubit case can be cast here 
in terms of appropriate Langevin equations for the quadrature 
operators of modes 1 and 2. Mode 1 is affected by white noise 
entering the cavity and mode 2 experiences quantum Brown- 
ian motion due to the contact of the mirror's vibrational mode 
with background phonons fl3ll . An analytical formalism for 
the treatment of this situation has been developed lfl3ll and 
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FIG. 4: (a): S(y) vs. v and n for the cbf model with J = 0.3, uj = 
0, n = 3, 7 = T = 0.001. (b): Same plot, for 7 = T = 10" 4 . The 
behavior against 7 is similar to Fig. [3](b). 



used to show that, by measuring the spectrum of the field, 
complete information regarding the temperature and damping 
rate of the mirror can be gathered 0, 0. For our purposes, 
it is important to stress that in the working conditions of the 
experiment in |4|], no entanglement is set between modes 1 
and 2, therefore putting the framework of Ref. [13] within the 
context of environmental-probing without entanglement. This 
case is doubly relevant as it describes a detecting stage (the 
field) also affected by environmental effects (the white noise 
entering the cavity), therefore demonstrating the validity of 
our approach for non-unitary dynamics of the "detector". 

Remarks.- In the characterization of the environment affect- 
ing a partially-accessible register, entanglement shared with 
a detecting device is not crucial. Even classical correlations, 
set through a time-independent interaction, help in identifying 
the nature of open- system dynamics and we have shown how 
an experimentally achievable figure of merit can be used. This 
result, achieved by studying the stationary state of the register- 
detecting stage system, is independent of the dimension of the 
Hilbert space and holds in a minimum-control scenario. We 
used our experimentally-oriented study to reinterpret the re- 
cent achievement of control in some mesoscopic systems |0] 
as genuine interaction-enabled environmental probing. 

Acknowledgments.- MP thanks C. Di Franco for discus- 
sions. We acknowledge support from the UK EPSRC, the 
QIPIRC and The Leverhulme Trust (ECF/40157). 



[1] A. Wallraff atal, Nature (London) 431, 162 (2004); D.I. Schus- 
ter et al, ibidem 445, 515 (2007). 

[2] RH.L. Koppens, etal, Nature (London) 442, 766 (2006). 

[3] M.D. LaHaye, et al, Science 304, 74 (2004); A. Naik, et al, 
Nature (London) 443, 193 (2006). 

[4] S. Gigan et al Nature (London) 444, 67 (2006). 

[5] L. Venkataraman, et al, Nature (London) 442, 905 (2006); J. 
Park, et al, ibidem 417, 722 (2002). 

[6] M. Greiner, etal, Nature (London) 415, 39 (2002). 

[7] I. Rau, G. Johansson, and A. Shnirman, Phys. Rev. B 70, 
054521 (2004). 

[8] See J. Lee, et al Phys. Rev. A 70, 024301 (2004); D. Ahn, 
J. Lee, M.S. Kim, and S.W. Hwang, Phys. Rev. A 66, 012302 
(2002). 

[9] A. Imamoglu, et al, Phys. Rev. Lett. 83, 4204 (1999); M. Pa- 
ternostro, et al, Phys. Rev. A 71, 042311 (2005). 

[10] M. Paternostro, G. Falci, M. Kim, and G.M. Palma, Phys. Rev. 
B 69, 214502 (2004); I. Wilson-Rae and A. Imamoglu, Phys. 
Rev. B 65, 235311 (2002). 

[11] D.F. Walls and G.J. Milburn, Quantum Optics (Springer, Hei- 
delberg, 1994). 

[12] J. Lee etal, J. Mod. Opt. 47, 2151 (2000). 

[13] M. Paternostro, etal, New J. Phys. 8, 107 (2006). 

[14] S.M. Barnett and P.M. Radmore, Methods in Theoretical Quan- 
tum Optics (Oxford, New York, 1997). 

[15] M.S. Kim and N. Imoto, Phys. Rev. A. 52, 2401 (1995); H. 
Jeong, J. Lee and M.S. Kim, Phys. Rev. A 61, 052101 (2000). 

[16] S. Bose, K. Jacobs, and PL. Knight, Phys. Rev. A 56, 4175 
(1997). 



